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ABSTRACT 


Target  component  feature  extraction  is  an  area  of  considerable  importance  to  the  Ballistic  Missile 
Defense  (BMD)  community.  In  particular,  extracting  essential  target  features  from  measurement  data  on 
targets  of  interest  leads  to  potential  target  identification.  The  extracted  component  features  correspond  to 
a  numerical  characterization  of  the  Geometrical  Theory  of  Diffraction  (GTD)  diffraction  coefficient  and 
can  also  be  used  to  develop  a  computationally  efficient,  measurement  based  RCS  signature  prediction 
model.  A  key  attribute  of  the  resulting  computational  model  is  that  the  measured  RCS  is  directly 
incorporated  into  the  computational  model. 

An  essential  ingredient  in  forming  a  measurements-based  signature  model  valid  over  a  wide  range 
of  frequencies  and  angles  is  the  ability  to  map  the  field  measurement  data  (2  dimensional)  onto  a 
component-based  three-dimensional  (3D)  geometry.  To  accomplish  this,  3D  characterization  of  the  target 
scattering  components  is  required.  Typically,  this  3D  characterization  of  the  scatterer  locations  is  obtained 
by  forming  a  3D  image  of  the  target,  and  extracting  the  dominant  scattering  centers.  In  this  paper  we 
extend  the  novel  formulation  for  3D  radar  imaging  of  Inverse  Synthetic  Aperture  Radar  (ISAR)  sparse- 
angle  data  using  high-resolution  spectral  estimation  theory  presented  in  a  previous  paper  (Ref  1)  to  the 
special  case  where  one  has  apriori  information  about  the  target  geometry. 

For  some  applications  one  has  apriori  information  about  the  target  geometry — for  example  in 
processing  static  range  data  or  in  processing  flight-test  field  data  on  ‘"known”  targets.  In  this  paper,  we 
develop  a  methodology  which  incorporates  this  known  surface  geometry  constraint  to  generate  a  3D 
ISAR  image  identifying  the  principal  scatterers  contributing  to  the  target  radar  cross-section.  The 
methodology  uses  only  a  sequence  of  single  snapshot  data,  from  which  key  features  on  the  target  can  be 
readily  extracted.  The  technique  eliminates  the  need  for  a  range-Doppler  correlator  required  to  track  the 
motion  of  any  specific  scatterer.  This  is  accomplished  by  hypothesizing  a  given  target  surface  constraint 
(3D),  and  mapping  sequential  2D  range-Doppler  location  estimates  onto  this  surface.  The  process  is 
adaptive,  as  the  assumed  geometry  is  modified  as  appropriate  (e.g.,  for  intelligence  related  apriori 
information)  to  match  the  data.  When  applied  to  the  sequence  of  high-resolution  snapshots  obtained  using 
the  direct  spectral  estimate  obtained  from  sequential  blocks  of  data,  it  presents  a  computationally  efficient 
way  of  generating  a  high-resolution  3D  image  from  which  the  dominant  scattering  centers  on  the  given 
surface  constraint  are  identified  and  characterized. 
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1.  INTRODUCTION 


Target  component  feature  extraction  is  an  area  of  considerable  importance  to  the  Ballistic  Missile 
Defense  (BMD)  community.  In  particular,  extracting  essential  target  features  from  measurement  data  on 
targets  of  interest  leads  to  potential  target  identification.  The  extracted  component  features  can  also  be 
used  to  develop  a  computationally  efficient,  measurement-based  RCS  signature  prediction  model.  A  key 
attribute  of  the  resulting  model  is  that  the  measured  RCS  is  directly  incorporated  into  the  computational 
model. 

An  essential  ingredient  in  forming  a  measurements-based  signature  model  valid  over  a  wide  range 
of  frequencies  and  angles  is  the  ability  to  map  the  field  measurement  data,  2  dimensional  (2D),  onto  a 
component-based  three-dimensional  (3D)  geometry.  To  accomplish  this,  3D  characterization  of  the  target 
scattering  components  is  required.  Typically,  this  3D  characterization  of  the  scatterer  locations  is  obtained 
by  forming  a  3D  image  of  the  target,  and  extracting  the  dominant  scattering  centers.  In  this  paper  we 
extend  the  novel  formulation  for  3D  radar  imaging  of  Inverse  Synthetic  Aperture  Radar  (ISAR)  data 
using  high-resolution  spectral  estimation  theory  presented  in  a  previous  paper  (Ref  1)  to  the  special  case 
where  one  has  some  apriori  information  about  the  target  geometry. 

In  previous  papers,  we  developed  a  novel  formulation  of  the  3D  imaging  problem  based  on  a  simple 
mathematical  representation  for  the  2D  to  3D  image  transformation,  which  processes  a  sequence  of 
independent  2D  high-resolution  snapshots  of  the  target  to  generate  the  3D  image  points.  The  formulation 
uses  modern  2D  spectral  estimation  techniques  to  enhance  the  resolution  obtained  using  conventional 
correlation  imaging  techniques.  This  new  technique  offers  several  advantages: 

1 .  Provides  for  a  direct  3D  estimate  (vs.  back  projection  to  a  3D  target  grid  matrix)  of  the  location 
of  the  dominant  scattering  centers,  having  enhanced  resolution  both  in  range  and  cross-range 
relative  to  conventional  imaging. 

2.  Provides  a  direct  estimate  of  each  scatterer  amplitude  and  phase  as  a  function  of  aspect-angle  to 
the  target. 

3.  Produces  a  3D  image  with  only  a  minimum  number  of  independent  snapshots,  where  each 
snapshot  represents  a  high-resolution  range-Doppler  image  at  a  known  aspect  angle.  Because  of 
the  snapshot  nature  of  the  technique,  it  is  particularly  applicable  to  3D  imaging  of  sectors  of 
sparse-angle  data,  for  which  the  sidelobes  of  the  correlation  imaging  integral  become  high. 

4.  Gracefully  reduces  to  a  range  only  (e.g.  wideband  radar,  low  pulse  repetition  frequency  (PRF)) 
and  Doppler  only  (e.g.  narrowband  radar,  high  PRF)  3D  imaging  capability. 

5.  By  implementing  a  range-Doppler  correlation  tracker  on  a  sequence  of  pulses,  one  can  isolate 
the  3D  motion  of  any  specific  scatterer  using  the  simple  2D  to  3D  mapping  developed  later  in 


the  paper,  and  extract  the  scatterer  amplitude  and  phase  as  a  function  of  time.  Using  this 
scatterer  response,  combined  with  the  3D  scatterer  locations  extracted  from  the  3D  image,  one 
can  directly  reconstruct  the  actual  scattering  measurements  using  a  simple  point  scatter  model 
based  on  the  Geometrical  Theory'  of  Diffraction  (GTD).  Correlating  the  extracted  GTD  based 
diffraction  coefficients  to  hypothesized  characteristics  of  the  scatterer  response  provides  the 
potential  of  enhanced  scatterer  ty  ping  and  identification. 

For  some  applications  one  has  apriori  information  about  the  target  geometry'  -  For  example  in 
processing  static  range  data  or  in  processing  flight-test  field  data  on  known  targets.  Of  prime  interest  in 
this  case  is  identification  of  the  dominant  scattering  centers  on  the  target,  with  the  intent  to  form  a 
computationally  efficient,  measurements  based,  signature  model  of  the  target.  In  this  report,  we  develop  a 
methodology  which  incorporates  this  known  surface  geometry  constraint  to  generate  a  3D  ISAR  image 
identifying  the  principal  scatterers  contributing  to  the  target  radar  cross-section.  The  methodology  uses 
only  a  sequence  of  single  snapshots  of  data  from  which  key  features  on  the  target  can  be  readily  extracted. 
The  technique  eliminates  the  need  for  a  range-Doppler  correlator  required  to  track  any  specific  scatterer. 
This  is  accomplished  by  hypothesizing  a  given  target  surface  constraint  (3D),  and  mapping  sequential  2D 
range-Doppler  location  estimates  onto  this  surface.  The  process  is  adaptive,  as  the  assumed  geometry  is 
modified  as  appropriate  to  fit  the  measured  data.  When  applied  to  the  sequence  of  high-resolution 
snapshots  obtained  using  the  direct  spectral  estimate  obtained  from  sequential  blocks  of  data,  it  presents  a 
computationally  efficient  way  of  generating  a  high-resolution  3D  image  where  the  dominant  scattering 
centers  on  the  given  surface  constraint  are  identified  and  characterized. 

Because  the  technique  provides  an  independent  estimate  of  the  location  of  each  scattering  center, 
and  a  separate  aspect  angle  and  frequency  dependence  for  this  scatterer,  when  we  refer  to  an  '"image'’  in 
this  paper,  we  consider  only  the  positions  of  the  scatterers  in  target  centered  space.  A  conventional  image 
also  provides  a  measure  of  the  average  intensity  at  these  locations,  usually  plotted  on  a  color  intensity 
scale  on  the  3D  grid  matrix.  However,  for  real  targets,  this  complex  amplitude  varies  considerably  with 
look  angle,  and  we  chose  to  delineate  this  variation  separately.  It  is  in  fact  this  component 
characterization  which  leads  to  the  measurements  based  signature  model  of  the  target. 
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2.  2-DIMENSIONAL  ALL-POLE  SIGNAL  MODELING  (REVIEW) 


A  general  development  of  the  3  D  signal  model  characterizing  the  high-frequency  field  scattered 
from  a  specific  target  using  the  Geometrical  Theory  of  Diffraction  (GTD)  and  its  relationship  to  2D,  all¬ 
pole  signal  modeling  techniques  was  developed  in  Reference  1.  For  completeness,  we  repeat  this 
development  in  this  section.  We  assume  a  target  centered  coordinate  system  as  illustrated  in  Figure  I . 


Figure  /.  Coordinate  System  and  Notation 


Standard  spherical  coordinates  0,  (j)  characterize  the  look-angle  to  the  radar,  which  is  determined  by 
the  unit  vector  k  given  by 


k  =  sin  ^cos  (/>k  sin  ^sin  (/>y  cos  0z  (1 ) 

where  the  notation  ‘V  denotes  a  ‘‘unit  vector.”  The  target  coordinates  are  defined  relative  to  the 
origin  of  this  system,  and  an  arbitrary  scattering  point  on  the  target  is  denoted  as  (x„,  yn,  z,i).  We  assume 
the  motion  is  known,  so  that  in  the  target  centered  reference  frame,  the  look  angle  of  the  radar  changes 
parametrically  with  time;  i.e.,  9  =  9  (t),  <])  =  (j)  (t). 
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Using  the  Geometrical  Theory  of  Diffraction  (GTD)  the  field  En  scattered  from  an  arbitrary  point 
located  at  (Xn,  yn,  z„)  is  given 


(2) 


where  c=  3  x  10^  m/s,  f  denotes  frequency,  denotes  the  vector  location  of  (Xp,  yp,  2^,)  and  D„ 
(/,  (fi)  is  the  well-known  scattering  diffraction  coefficient  introduced  by  Keller  [Ref  2].  In  a 
generalized  sense,  the  scattering  coefficient  characterizes  a  broad  class  of  scattering  features.  It  will  be 
seen  later  that  it  corresponds  directly  to  the  amplitude  and  phase  extracted  from  measurement  data  using 
2D  all-pole  spectral  estimation  techniques. 

In  general,  D„  by  itself  is  a  complex  function  whose  frequency  and  angle  variation  is  considerably 
slower  than  the  total  phase  variation  in  (2).  Thus  a  snapshot  of  data  provides  information  only  on  the  local 
variation  of  D„  when  viewed  as  a  wide  angle,  broadband  function. 

Wide-angle  characterization  of  Dp  is  accomplished  by  processing  a  '‘block”  of  pulses  consistent 
with  changes  in  0  (t),  ^  (t).  Eq  (2)  is  rigorously  valid  in  the  frequency  domain,  for  a  single  pulse.  For  the 
time  scales  considered  here,  the  time  variable  “t”  represents  “pulse-to-pulse'’  changes  in  the  radar  return. 
Since  the  changes  in  0,  (j)  are  not  directly  radar  observables,  the  radar  return  pulses  are  processed  as  a 
function  of  frequency  and  time.  In  frequency-time  space,  the  measurement  of  E„  can  be  expressed  in  the 
form: 


(3) 


where  vj/p  (f,  t)  and  Dp  (f,  t)  are  implicitly  a  function  of  k  and  according  to 


4  7T  f  ^ 

Vn  (f’  0  =  -  k(t)-  r„ 

c 

D,(f.  t)  =  D„(f,  k-  rj 


(4) 


Equations  (3)  and  (4)  provide  projections  defined  relative  to  a  3D  signal  model.  However,  the  radar 


observables  (f,  t)  comprise  a  2D  space.  Thus  a  mapping  from  2D  radar  observable  space  to  3D  target 
centered  space  (i.e.;  Xp,  yp,  Zp)  is  required  for  3D  imaging. 
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Consider  the  local  approximations  of  y/ ^  and  Dn  over  frequency  and  time.  From  the  concept  of  all¬ 
pole  modeling,  Dn  (f,t)  can  be  locally  approximated  as  having  an  exponential  behavior,  i.e.. 


D,.  (f,  t)~Do„  e- 


„f  e-P„. 


(5) 


where  an  and  Pn  characterize  the  localized  exponential  approximation  and  Don  is  a  complex 
constant.  This  local  approximation  is  particularly  convenient  for  using  2D  all-pole  modeling  spectral- 
estimation  techniques,  for  which  the  poles  characterizing  the  signal  model  have  an  exponential  behavior 
with  complex  exponent.  The  phase  behavior  (f,  t)  is  somewhat  more  complicated.  Consider  two 
special  cases: 


a.  Fixed  0^^  (j>^.  Vary  f  =  fo  +  kAf?  k  =  l,...K 
Then  we  can  expand  y/ ^  (f,  to)  in  the  form 


Wn  (f’  t0)=  Wxv  (fO’  to)  + 


5f 


A/ 


where  to  corresponds  to  Go  =  0  (to),  (j)o,  =  <|)  (to). 


(6) 


b.  Fixed Vary  t  =  tQ  +  q  A  t,  q  =  1,...,  Q 

For  this  case  we  can  expand  ij/ ^  (fo,  t)  in  the  form 


0~t^n  (fo’  to)  + 


di 


q 


(7) 


If  Af  and  At  are  small  enough,  Eqs  (6)  and  (7)  can  be  combined  to  form  a  2D  all-pole  signal  model 
valid  for  linear  imaging  (e.g.,  as  for  the  conventional  range-Doppler  image).  In  this  case,  the  field  E„  (f,  t) 
becomes  a  function  of  the  discrete  2D  sequence,  En  (k,  q),  which  can  be  written  in  the  form 
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(8) 


E„  (k,  q)  =  p- 


where  the  complex  poles  s„  and  pn  are  given  by 


(9) 


Pn  =  e 


P„+J 


(10) 


A  conventional  linear  image  could  then  be  obtained  using  a  weighted  2D  Fourier  transform  on  Fq 
(8).  However  enhanced  resolution  is  realized  when  the  poles  Sn,  Pn  are  determined  from  a  2D  pole 
estimation  algorithm  applied  to  the  measurement  data  set.  Once  the  pole  pairs  (s,,,  Pn)  are  determined,  the 
complex  amplitude  Do,,  can  be  extracted  using  a  "least  squares"  fit  of  a  sum  of  poles,  each  weighted  by 
Don  in  the  form  of  Eq  (8),  to  the  measured  data  set. 

Figure  2  illustrates  the  methodology  used  for  sequential  processing  of  a  block  of  uncompressed 
radar  data  pulses.  The  poles  s„  and  pn  are  extracted  from  the  composite  scattered  field  data  matrix  E  (k,  q) 
and  converted  to  range  and  range  rate  space  using  (9)  and  (10).  We  obtain 


(II) 


Using  Eq  (4),  we  obtain 


dy/„ 


df 


(12) 


r  n  —  2— 
dt  c 
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where  we  define  k  =  dk  /  dt.  Substituting  Eq  (12)  into  (1 1),  we  obtain  the  desired  mapping: 


K  =  k-  r„  = 


4;rAf 


Im  (InsJ 


R„  =  k  •  r„  = 


4;r  A  t 


Im  (In  pj 


(13) 

(14) 


where  we  define  R^)  to  be  the  range,  range  rate  observables  corresponding  to  (Xn,  yn,  Zn)  and 
Im  (z)  denotes  the  ‘‘imaginary  part”  of  the  complex  variable  z.  The  complex  amplitude  associated  with 
each  spectral  component  (scatterer)  is  computed  by  the  2  D  spectral  estimation  algorithm.  The  sequence 
{R„,  R„,  Dn}q  can  then  be  used  for  image  interpretation  and  component  extraction  and  characterization. 
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3.  SEQUENTIAL  BLOCK  PROCESSING  FOR  COMPONENT  EXTRACTION 

AND  CHARACTERIZATION  (REVIEW) 


The  sequence 


{Rn.  Rn.  D„}, 


where  n  denotes  the  scatterer  present  in  the  data  at  time 


snapshot  t  =  tq,  provides  a  sequential  characterization  of  the  target  scattering  components  and  RCS  as  a 
function  of  time.  We  now  show  how  this  sequence  can  be  used  for  image  interpretation  and  component 
extraction  and  characterization. 


Component 

/ 

. :z 

f 

I 

Motion 

Image 

Extraction 

! 

Estimation  ^ 

[  Interpretation  , 

and 

Characterization  j 

Range,  Range  Rate  Sequence 
For  each  resolvable 
Scatterer:  n  =  1, N 


Amplitude,  Phase 
Sequence  for  each 
Resolvable  scatterer 


Figure  2 


Figure  2.  Sequential  Block  Processing 


Figure  3  illustrates  a  typical  result  for  a  single  data  block  using  simulated  data  on  a  spinning 
cylinder  with  ‘‘fins”  attached  at  one  end.  Two  specific  2D  spectral  estimation  algorithms  have  been 
developed,  one  based  on  ESPRIT  processing  techniques  (Ref  3)  and  the  other  based  on  a  state-space 
formulation  (Refs  4,5).  The  resulting  (R„,  R  „)  estimates  are  overlaid  on  a  conventional  weighted  Fourier 
Transform  based  range-Doppler  image.  The  ®  marked  points  on  the  image  represent  the  direct  estimate  of 
the  scatterer  locations  to  within  the  accuracy  of  the  Cramer-Rao  bound  obtained  by  applying 
ESPRIT  directly  to  the  data  block.  Processing  these  high-resolution  (/?^,  estimates  provide 
several  advantages  over  processing  the  conventional  image: 
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1 .  The  estimate  of  (  /?,  R  )  for  each  scatterer  is  direct  -  i.e.  it  need  not  be  inferred  from  the  image. 
The  resultant  extraction  provides  a  computationally  efficient  numerical  sequence  for  further 
processing  to  determine  motion  and  component  characterization. 

2.  The  resolution  and  accuracy  achieved  closely  approaches  the  Cramer-Rao  bound  for  each 
scatterer  location. 

3.  The  amplitude  and  phase  of  each  scatterer  is  a  direct  output  of  the  spectral  estimate,  and  need 
not  be  inferred  from  the  image  (such  as  when  using  the  so-called  ‘‘clean”  extraction  algorithm 
applied  to  a  conventional  DFT-based  range-Doppler  image).  For  example,  the  3D  image 
generated  can  be  constructed  using  only  scatterers  above  a  certain  amplitude,  thus  eliminating 
“noisy”  or  unimportant  scatterers. 


Data  Block 


time 


Conventional 

FFT 

Weighted 

Image 


2  D  Block 
Spectral  Estimate 

(R„,  R„,  D„) 

<8) 


Simulated>Cylinder  with  Fins 


-ISO  -100  -so  □  S]  !□□  ISO 


Figure  3 


- ► 

Range  Rate 


Figure  3.  High-resolution  Block  Processing  v^.  Conventional  DFT  Image 
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4.  HIGH-RESOLUTION  3D  IMAGING  USING  CUMULATIVE  SNAPSHOTS 


The  technique  developed  in  Reference  [I]  using  cumulative  snapshots  to  form  a  3D  image  will  be 
briefly  summarized  in  this  section.  It  is  most  useful  for  simpler  targets  characterized  by  a  small  number 
( £  10)  of  scatterers).  Each  pulse  return  is  uncompressed  and  de-weighted  to  obtain  the  frequency  domain 

representation  of  the  scattered  field  for  that  particular  look  angle.  A  2D  spectral  estimation  algorithm  is 
applied  to  a  block  sequence  of  pulses  to  estimate  Sn,  Pn.  and  Don.  A  /?,  R  correlator/smoother  is 
implemented  to  track  these  poles  in  /?,  R  space,  sequentially  sliding  (at  fixed  jump  intervals)  a  block  of 
pulses  through  the  measurement  data  set.  Eq  (13)  and  (14)  are  then  applied  sequentially  to  a  number  of 
snapshots  to  develop  an  estimate  of  the  scattering  points  (Xn,  yn,  Zn)-  To  illustrate  this  process 
mathematically,  we  express  the  transformation  defined  by  (13)  and  (14)  in  matrix  form: 


T  • 

=^i 


Xn 

yn 

— 

K. 

-^1  . 

(15) 


where  T  is  a  2  x  3  transformation  matrix  given  by 


T  = 


(16) 


where  k  and 
readily  shown  that 


k  are  expressed  as  row  vectors,  and  Eq  (15)  defines  the  mapping  at  t  —  tci\.  It  is 


k  - 


(17) 


We  can  accumulate  these  equations  for  two  snapshots  at  tqi  and  tq2  in  the  form 


(18) 


[T  1 

(x  ^ 

=<1, 

yn 

L 

^2  _ 

\^n  ) 

\^n  ) 

(K'\ 


\^nj 


^2 


As  additional  snapshots  are  obtained,  these  equations  are  added  to  (18),  which  is  then  solved  in  a 
least  squares  sense  for  the  scattering  coordinates  (x„,  Zn). 

Note  that  the  accumulation  of  snapshots  emphasizes  the  importance  of  a  good  /?,  R  tracker, 
particularly  in  low  signal-to-noise  ratio  environments.  It  is  imperative  in  Eq  (18)  that  the  sequence  of 
Rp  estimates  correspond  to  the  same  physical  scatterer  located  at  (Xn,  yn,  z,,).  For  the  results 
presented  in  Ref  1,  we  used  a  simple  "nearest  neighbor"  tracking  algorithm  in  R,  Rdot  space  to  track  the 
motion  of  the  scatterers. 

In  the  next  section,  we  develop  a  simple  version  of  (18)  not  requiring  an  accumulation  of  snapshots, 
by  adding  a  third  equation,  obtained  by  imposing  a  surface  constraint  to  the  two  equations  defined  by  Eq 
(15). 
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5.  3D  SNAPSHOT  IMAGING  USING  A  SURFACE  CONSTRAINT  - 

DEVELOPMENT 


Expand  the  transformation  matrix  T  into  a  more  general  form  defined  by 


T  T  T 

^11  ^12  ^13 

T  T  T 

/2l  ^22  ^23 


(19) 


Thus,  each  scatterer  (Xn,  yn,  z„)  on  the  surface  S  must  satisfy  the  condition 


Til  x„  +  Ti2y„  +  Ti3Zn  =  R„ 


(20) 


T21  x„  +  T22  yn  +  T23  z„  -  (21) 

If  the  surface  S  is  known  (or  hypothesized),  then  an  additional  equation  completing  (20)  and  (21) 
above  can  be  obtained  by  requiring 


S  (Xn,  yn,  Zn)  =  0 


(22) 


where,  in  general,  S  (x,  y,  z)  =  0  describes  the  target  surface  geometry.  Imposing  the  condition  Eq 
(22)  on  (20)  and  (21)  provides  a  direct  set  of  generally  independent  equations  which  define  a  (x„,  y,,,  Zn) 
scatterer  position  using  a  single  snapshot  of  data.  The  tradeoff  that  we  will  elaborate  on  later  is  the 
existence  of  several  types  of  ambiguous  solutions,  which  map  to  the  surface  but  do  not  correspond  to 
physical  scatterers  on  the  target.  We  address  this  topic  further  in  section  6. 

Two  simple  cases  will  be  considered  in  this  section  for  example  purposes:  a  cone  and  a  cylinder. 
These  provide  for  a  simple,  direct  solution  for  constructing  a  3D  image  directly  from  the  field 
measurement,  which  can  be  “thresholded"’  by  prefiltering  on  the  amplitudes  |Dn|  associated  with  each 
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scatterer.  Although  simplistic  in  nature,  the  examples  lead  to  a  quantitative  understanding  to  the 
application  of  the  technique.  For  the  two  cases  considered,  we  have 

a.  Cylinder  (Radius  a) 


S  (x,  y,  z)  =  x^  +  y^  -  =  0  (23) 

b.  Cone  (Cone  half-angle  0^ ) 


S  (x,  y,  z)  =  z -  Zj,  +  ctn  0^  •  -Jx^  +y^  =  0  (24) 

where  ‘V'  is  the  cylinder  radius  and  Zb  denotes  a  surface  bias  (shift)  used  to  align  the  cone  with  the 

data. 

Eqs  (20)  —  (22)  although  nonlinear  in  nature,  can  be  solved  exactly  for  these  two  simple  surface 
constraints.  For  both  cases,  Eqs  (20)  and  (21)  can  be  solved  to  express  x„  and  y„  as  functions  of  Zn  and 
written  in  the  form 


Yn  =a,  z„  +  b, 


(25) 


z„  +  h, 


(26) 


a,  = 


T  T 

^13  ^2! 


-  T  T 

^11  ^23 


(27) 
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(28) 


b.  =  - 


T,,  R,  -  T„  R, 


a,  =  - 


T,,  Tj;  -  T,;  T 
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R,  T„  -  T,,  R, 


(29) 


(30) 


and 

A  —  T,i  T22  -  T,2  T2J 


(31) 


Using  (25)  and  (26)  for  Xn  and  y„  in  either  (23)  or  (24)  defining  S  (x,  y,  z)  produces  a  quadratic 
equation  for  Once  Zn  is  determined,  Xn  and  yn  can  be  computed  from  (25)  and  (26).  Note  that  since  (23) 
and  (24)  are  quadratic  in  form,  two  solutions  for  Zn  will  be  generated.  For  the  simple  cases  considered 
here,  the  correct  solution  is  readily  determined.  However  in  the  more  general  case,  these  ambiguous 
solutions  must  be  filtered  out  as  discussed  in  Section  6. 

5.1  A  SIMPLE  EXAMPLE  (SIMULATION) 

To  illustrate  the  technique,  consider  the  simulation  considered  in  Ref  [1],  and  illustrated  in  Figure  4. 
The  target  consists  of  a  spherical  tipped,  conical  model  of  a  re-entry  vehicle  (RV)  of  length  =  1.6  m, 
having  4  antenna  ports  located  at  90°-roll  intervals  and  having  a  flat  base.  Each  antenna  port  was 
simulated  to  be  a  constant  amplitude  scatterer  independent  of  roll  and  5  dB  reduced  from  the  nose.  We 
simulated  wideband  pulses  (using  standard  Physical  Theory  of  Diffraction  (PTD)  codes  to  compute  the 
net  scatter  vs.  frequency,  and  a  hamming  weighted  Fourier  Transform  over  frequency  to  generate  the 
compressed  pulse  shapes)  received  by  an  X-band  sensor  for  motion  defined  by  a  spinning  target  without 
coning.  Figure  4  illustrates  the  RV  geometry  and  target  orientation  for  0o  =  40°.  An  X-band  bandwidth  of 
1  GHz  is  used  for  example  purposes.  Data  is  processed  for  a  variety  of  snapshots  taken  over  a  360°  roll 
cycle.  Also  shown  in  Figure  4  are  the  composite,  weighted  pulse  shapes  received  at  the  sensor  over  the 
full  360°  roll  cycle,  where  the  radar  line  of  sight  to  the  target  is  from  right  to  left  —  i.e.,  the  nose  is  roughly 
10  dB  larger  than  the  base.  Note  that  the  antenna  ports  are  clearly  resolvable  from  the  nose  and  base.  For 
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any  given  roll  angle,  two  antenna  ports  are  visible  in  the  same  range  resolution  cell,  resulting  in  a  pulse- 
pulse  variation  of  the  pulse  return  in  the  RV  midsection  due  to  the  interaction  of  the  two  antenna  port 
scatterers.  Since  the  nose  and  base  of  the  target  are  roll  symmetric,  the  pulse  returns  from  these  scatterers 
are  insensitive  to  time  for  constant  aspect  angle  and  hence  these  scatterers  produce  zero  Doppler  return. 
Thus  they  correspond  to  so-called  ‘"slipping  scatterers.” 


RV  with  4  Antenna  Ports 
(includes  shadowing) 


X-Band  Composite  Pulse  Shapes 
(I  Ghz  l^andwidth) 


Figure  4 


Figure  4.  Constant  Aspect  Angle  Imaging  (Spinning  RV  with  Shadowed  Antenna  Ports) 


Figure  5  illustrates  the  {R^,  R^}  sequence  extracted  from  the  data  over  one  complete  roll  cycle 


of  the  target,  and  the  resulting  3D  RV  “image”  reproduced  from  Ref  1  using  the  (Xn,  yn,  z^)  computed  by 
correlating  the  scatterers  and  accumulating  snapshots  according  to  Eq  (18).  Note  that  a  multi-aspect  wide- 
angle  image  is  required  for  complete  target  imaging  because  of  shadowing.  For  the  pure  spinning  motion 
considered,  the  base  of  the  target  presents  zero  Doppler  to  the  radar  line  of  sight,  and  thus  appears  as  a 
slipping  scatterer.  For  the  image  illustrated,  we  modify  the  conventional  image  of  the  base  return  as  a 
slipping  scatterer  by  moving  the  base  slipping  scatterer  to  its  projected  (x„,  yn,  Zn)  location  assuming  an  a 
priori  known  base  diameter.  This  model-based  approach  to  imaging  slipping  scatterers  is  often  used  when 
additional  data  on  the  target  from  known  sources  or  other  wide-angle  aspect  angle  observations  is 
available.  Alternately,  if  motion  change  in  the  orthogonal  direction  is  available,  e.g.,  if  the  RV  is  “coning’’ 
as  well  as  spinning,  additional  orthogonal  snapshots  could  be  added  to  the  image  processing. 


Figure  6  illustrates  the  3D  image  generated  using  the  conical  surface  constraint  defined  by  Eq  (24) 
for  0^  =  7.3  .  Gaps  in  the  slipping  scatterer  image  appear  because  the  noise  inherent  in  the  R  poles 
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that  occur  as  each  antenna  port  approaches  the  shadow  region  were  filtered  out.  Note  also  that  the  base 
slipping  scatterer  is  automatically  projected  to  the  proper  position  on  the  cone  surface,  and  is  consistent 
with  the  model-based  projection  used  in  Figure  5  assuming  an  apriori  known  base  diameter. 

Figure  7  illustrates  the  sensitivity  of  the  images  to  cone  angle  0^^  as  indicative  of  using  the  cone 
angle  as  a  parameter  to  "‘focus”  the  image  when  the  cone-angle  might  not  be  apriori-knovvn  exactly. 
Clearly  best  results  occur  for  0^  =  7.3  ,  where  the  resulting  variance  in  antenna  port  locations  is  due  to 
the  variance  in  estimate  shown  in  Figure  5. 


{RJq  Sequence 

3D  -  RV  “Image” 


►  time 
Figure  5 


Figure  5.  Cumulative  Block  Processing  to  SD-Target  "Image"' 
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Antenna  Ports 


t  O  f 

v-y 


meters  Figure  6  meters 

Figure  6.  3D  Snapshot  Image  with  Sequential  Conical  Surface  Constraint 


Antenna  Ports 

C  e.  =  8° 

I  0  1  ^ 

v-y 

Top  View 
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meters  Figure  7  meters 

Figure  7,  Image  Defociising  vs.  Error  in  Cone  Angle  (Top  View) 
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6,  A  GENERAL  “CAD-MODEL”  SURFACE  CONSTRAINT 


The  examples  of  the  previous  section  demonstrate  the  utility  of  applying  the  surface  constraint  to  a 
simple  conical  RV-like  surface  geometry.  However  the  cone  and  cylindrical  surfaces  are  usually  too 
simplistic  to  be  useful  for  application  to  complex  targets,  and  a  more  general  approach  is  required.  We 
now  develop  this  more  general  approach  by  considering  the  surface  S  (x,y,z)  =0  to  be  defined  numerically 
at  a  fixed  number  of  grid  points 


S  (Xp,  yp,  Zp)  =  0,  p  =  I ,.  .  .,Ng  (32) 

where  Ng  denotes  the  total  number  of  grid  points  defining  the  target  surface.  Consider,  for  example, 
a  simplistic  model  of  a  missile  rocket-body  consisting  of  a  finite  cylinder  having  four  equally  spaced  fins 
at  the  rear.  A  point-wise  numerical  representation  of  this  surface  is  illustrated  in  Figure  8,  showing  a  3D, 
top  and  side  view.  (Note  the  x,  y  scales  in  Fig.  8  have  been  enlarged  for  clarity).  Because  of  the 
similarities  of  this  type  of  surface  characterization  to  a  computer-aided-design  (CAD)  model  of  the  target, 
we  adopt  that  nomenclature,  recognizing  that  as  defined  here  the  model  might  not  be  suitable  for  input  to 
other  numerical  electromagnetic  computation  codes.  In  fact  one  useful  interpretation  of  the  point-wise 
characterization  of  the  surface  useful  in  solving  Eqs  (20),  (21)  and  (32)  numerically  is  to  view  each  point 
as  a  potential  scatterer.  The  required  density  of  the  surface  grid  is  then  determined  by  the  resolution  of  the 
sensor  (in  range  and  cross-range,  given  by,  respectively,  c/FBW  and  >-/FBW,  where  FBW  is  the  fractional 
bandwidth  of  the  sensor  -  i.e.,  FBW  =  BW/fo,  where  BW  is  the  actual  bandwidth  (Hz)  and  fo  the 
operating  frequency). 


Figure  8.  Rocket  Body^  CAD  Model  for  Simulation 
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The  technique  we  use  to  solve  Eqs  (20,  (21)  and  (32)  is  ‘"hypothesis-test”  based.  For  a  given 
incidence  angle  (0,  (|)),  the  scattering  from  the  target  is  characterized  by  N  pole  pairs  (converted  to  range, 
range  rate  using  Eqs  (13)  and  (14)),  Rn}’  n  =  1,  .  .  .,N.  In  effect  we  will  test  each  pair 

(Rj^,  Rj^)  corresponding  to  a  given  angle  to  check  if  it  emanates  from  a  particular  point  on  the  grid 
surface.  The  result  usually  corresponds  to  a  physical  scattering  center  and  several  ambiguous  solutions 
(i.e.,  the  point  on  the  body  is  consistent  with  Eqs  (20),  (21)  and  (32),  but  is  not  one  of  the  scattering 
centers  contributing  to  the  scattered  field).  We  then  apply  a  series  of  “ambiguity  filters”  to  this  solution, 
leading  to  the  desired  solution  set  {x„,  yn,  z„}.  This  process  is  repeated  for  each  incidence  angle,  and  the 
resultant  point-set  corresponds  to  the  dominant  scattering  centers  on  the  target.  When  superimposed  on 
the  CAD  model,  they  appear  as  an  “image”  of  point-scatterer  locations. 

Quantitively,  the  process  is  described  in  the  steps  below: 

4.  Since  most  ballistic  targets  are  nearly  axially  symmetric,  it  is  useful  to  sort  the  grid  matrix  {Xj,, 
yg,  Zg}  in  increasing  Zg;  i.e.  we  assume 


Zg  =  sort  [Zg] 


(33) 


where  Zg  (1)  <  Zg...<  Zg  (Ng),  and  relabel  Xg  and  yg  accordingly.  This  allows  one  to  identify 

the  front  or  rear  of  the  target  by  its  location  in  the  ordered  Zg  string. 

5.  Consider  now  the  set  {Rn .  Rn  }  characterizing  the  scattering  at  a  specific  angle  (0,  (j)).  Assume 
that  the  scattering  associated  with  each  pair  emanates  from  Zgp,  where  the  subscript  p  denotes  p‘^^ 
grid  point.  If  this  hypothesis  were  true,  then  the  corresponding  (x,  y,  z)  scatterer  coordinates 
would  be  given  by  (x^^,  y^^,  Zgp)  ,  where  Xn  and  yn  are  determined  from 


(x  \ 


^  1  _  I 

^xy  -^gp  ^\y 


Rr 

V^ny 


(34) 


and  Txy  and  T^  are  defined  as 


'T,|T,2^ 

;  T,  = 

^^21  T22; 

’  z 

(35) 
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To  check  the  hypothesis,  we  compare  the  point  (xp,  Zgp  j  to  see  if  it  actually  lies  on  the  grid 
surface,  according  to  an  error  tolerance  Cj: 


|(xg-x„)^+(yg-yn)^|''^<sj 


(36) 


If  (36)  is  satisfied,  we  admit  the  point  (Xp,  yp,  Zgp)  as  a  possible  scattering  center. 

6.  Repeat  Step  2  for  each  grid  point  Zgp,  p  =  l,...Ng  and  store  the  collection  of  possible 
scattering  centers . 

7.  Apply  a  series  of  ambiguity  filters  to  the  initial  data  set  {xp,  yp,  Zp}.  These  filters  are 
developed  in  the  following  section. 

6.1  SCATTERING  MECHANISMS  AND  AMBIGUITY  FILTERS 

There  are  several  physical  mechanisms  for  which  ambiguous  scattering  centers  can  occur  using  the 
above  methodology.  Some  of  these  are: 

a.  “Shadowing” 

b.  Component  to  Component  Coupling 

c.  Slipping  Scatterers 

d.  CAD  Model  Surface  Tolerances 

We  now  develop  appropriate  “Ambiguity  Filters”  to  eliminate  these  ambiguous  solutions.  We  note 
that  the  specific  filters  to  use  are  strongly  target  (surface)  dependent,  and  those  discussed  here  are  focused 
on  the  cylinder-fin  surface  example  considered.  The  mechanisms  however  are  generally  applicable  to  any 
specific  target. 

6.1.1  Shadowing 

For  the  simple  cone  and  cylinder  examples  considered  in  Section  5,  multiple  roots  inherent  to  the 
quadratic  analytic  solution  developed  in  Section  5  lead  to  a  scatterer  located  in  the  shadow  region.  In  fact, 
although  not  developed  here,  it  is  possible  to  find  closed  form  expressions  for  (Xn,  yn,  Zn)  for  the 
cylindrical  constraint  (23)  for  the  case  0  =  0,  and  demonstrate  the  presence  of  one  of  the  quadratic  roots 
occurring  in  the  shadowed  region.  The  same  effect  will  be  true  for  the  more  general  case  of  the  cylinder- 
fin  CAD  example. 
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A  simple  check  can  be  imposed  to  filter  out  this  ambiguous  shadowed  solution.  Denote  the 
candidate  scatterer  point  as  (Xn,  yn,  Zn)-  The  scatterer  lies  in  the  cylinder  shadowed  region  whenever 


k-(Xn<yn’  0)<0 


(37) 


Eq  (37)  is  strictly  true  for  an  infinite  cylinder.  For  a  finite  cylinder,  scatterers  on  the  visible 
truncated  top  or  bottom  cylinder  surface  are  not  shadowed.  Thus  we  apply  (37)  only  to  grid  points 
Zg  <  z^> ,  where  Zc  is  the  height  coordinate  of  the  visible  cylindrical  cap. 

6.1.2  Component  to  Component  Coupling 

For  targets  that  have  some  components  having  large  radial  dimensions  (e.g.,  fins),  scattering  from 
the  outer  surface  can  lead  to  an  ambiguous  scatterer  projected  to  a  smaller  radial  diameter  surface  (e.g., 
cylinder).  Figure  9  illustrates  this  effect  for  the  fin-cylinder  surface  geometry.  (Note,  the  fins  are  only 
sparsely  populated  with  grid  points  in  this  example).  A  side-view  and  top  view  clearly  shows  how 
scatterers  having  the  same  (R,  R)  pairs  resulting  from  fin  scattering  can  project  to  the  smaller  diameter 
cylindrical  surface.  Thus,  a  fin-cylinder  coupling  ambiguity  filter  is  required  (Note,  to  emphasize  this 
effect  in  Figure  9,  we  have  chosen  the  error  tolerance  parameter  Ej  in  Eq  (36)  abnormally  large). 

It  is  easy  to  define  the  mathematical  condition  for  this  effect  to  occur.  Denote  (X|,  y|,  Z|)to  be 
the  true  fin  scatterer  location  and  (X2,  y2^  ’^2)  ambiguous  scatterer  location.  These  two  points 
satisfy  the  common  constraint 


^xy 


\y\J 


+  Z]  -  Tvv 


xy 


/'y  ^ 


+  Tz  Z2  - 


(38) 
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Error  tolerance  =  5cm 
(theta,phi)=(56J52)  deg 


Side  View 


Radar  LOS 


Figure  9.  Ambiguous  Scatterer  Location 


For  the  filter  design,  apply  (38)  to  the  condition  that  (xi,yi,zi,)  is  on  a  fin.  When  this  occurs, 
positions  (X2,  y2,  Zg)  where  Zg  is  an  arbitrary  grid  position,  must  satisfy 


X2 

vy2y 

+  T;J  (Z|  -  Zg) 


(39) 


We  simply  apply  (39)  for  every  grid  point  Zg  and  filter  out  those  points  lying  on  the  grid  cylindrical 
surface. 

6.1.3  Slipping  Scatterers 

Slipping  scatterers  are  unique  in  the  sense  that  they  exhibit  no  Doppler  with  change  in  roll  angle 
((j)).  Thus,  for  0  =  0,  R  =  0  for  a  slipping  scatterer.  Hence  by  definition  it  is  an  axially  symmetric 
scatterer  over  the  region  for  which  (j)  =  0 .  The  ambiguities  generated  for  a  slipping  scatterer  are  typically 
in  the  shadow-region,  and  are  filtered  out  by  this  filter.  However,  for  the  finite  cylinder,  the  back-side 
scattering  from  the  truncated  top  of  the  cylinder  (shadowed  for  the  infinite  cylinder)  is  clearly  visible  to 
the  radar  when  the  aspect  angle  is  less  than  90°.  The  ambiguity  associated  with  this  ‘'back-side’'  slipping 
scatterer  is  in  the  visible  region,  and  still  another  filter  is  required.  For  this  case,  since  it  is  unique,  we 
simply  compute  its  location  and  eliminate  the  point.  To  compute  its  position,  we  apply  (38)  assuming 
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(X|,  Yj,  Z|)  is  the  back-side  slipping  location  and  (X2,  y2,  Z2)is  the  visible  region  ambiguity.  To 
compute  Z2,  we  set  (j)  =  0,  (by  symmetry,  for  a  slipping  scatterer)  for  which  Xi  =  -a,  yi  =  0  and  z\  =  Zc, 
where  “a’'  is  the  cylinder  radius  and  is  defined  previously.  For  the  ambiguity,  X2  =  +  a,  y^  =  0.  Then  zi 
is  given  by 


Z2  =  z^.  -2a  tan9 


(40) 


The  X2,  y2  locations  are  the  mirror  images  of  the  back-side  slipping  scatterer,  i.e. 


X2  =  a  cos  (j),  y2  =  a  sin  ^ 


(41) 


For  the  filter,  we  eliminate  this  point  for  each  (0,  (j))when  the  cylinder  back  side  '"cap”  is 
illuminated. 

6.1.4  CAD  Model  Surface  Tolerances 

We  note  in  the  methodology  presented,  the  “surface  tolerance”  parameter  £j  defined  by  Eq  (36)  is 
key  to  determining  when  a  hypothesized  scatterer  lies  on  the  surface.  Typically,  if  £j  is  chosen  too  large, 
multiple  ambiguities  occur.  This  is  illustrated  in  Figure  9,  where  we  have  deliberately  chosen  Cj  large. 
Ideally,  Cy  is  chosen  on  the  order  of  the  grid-spacing  size.  Clearly,  there  is  a  trade-off  here.  If  Cj  is  chosen 
too  small,  and  the  grid  spacing  is  large,  no  scattering  centers  will  be  detected.  For  our  purposes,  Cj  is 
chosen  somewhat  empirically,  for  the  surface  geometry  considered. 
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7.  CYLINDER  -  FIN  TARGET  SIMULATION 


In  this  section,  we  simulate  the  back-scattered  field  from  a  cylinder-fin  type  ballistic  target  using 
standard  Physical  Theory  of  Diffraction  (PTD)  type  codes,  and  develop  a  CAD  numerical  surface  model 
to  illustrate  the  methodology  developed  in  Section  6.  We  assume  sensor  parameters  identical  to  those  for 
the  simulation  in  Section  5.1.  Figure  10  illustrates  the  missile  body  cylinder-fin  geometry  used,  as  well  as 
the  dynamical  motion  parameters  used  to  characterize  a  spinning  (fs=  spin  frequency),  precessing  (fp  = 
precession  frequency)  motion,  where  fs  »  fp  and  Tp  =  1/fp,  Ts=  1/fs.  Figure  1 1  illustrates  the  example 
motion  converted  to  (0,  coordinates,  and  a  range-time  intensity  plot  over  roughly  2  precession  cycles. 


The  sequence  |Rn,  Rn,  ^  ^  extracted  using  sliding  data  blocks  as  described  in  Section  3. 

The  resulting  sequence  for  (R^,  R^)  over  the  time  interval  [0.3,  0.5]  seconds  is  illustrated  in  Figure  12, 
and  the  sequence  for  |D„|  in  Figure  13.  Over  this  particular  time  interval,  the  radar  illuminates  the  top 
portion  of  the  cylinder.  As  can  be  seen  from  Figure  1 1,  the  aspect  angle  to  the  target,  over  this  time 
interval,  gradually  decreases  from  ~100  degrees,  through  90®  (broadside)  and  up  to  the  60®  minimum 
aspect  angle  determined  by  the  assumed  precession  motion.  Thus  the  observed  range  extension  of  the 
target  increases  as  a  function  of  time  (over  this  interval)  to  a  maximum,  and  then  begins  to  decrease  as  0 
increases  out  of  minimum  precession  angle. 


Motion  Parameters 

/c  =  120® 
ap  =  +75® 

as  =  0° 

Op  =  64® 

7>  = -0.455  sec 
7’j  =  0.0709  sec 
dt  =  4.05102^-5sec 


Figure  10.  Rocket-Body — CAD  Model  Simulation 
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Profile 


Range  (m} 


Simulated  Motion  (Spin  and  Precession) 


Figure  II.  RTI  Plot  of  Simulated  Rocket  Bod}' 


Order=10 


Order=5 


Figure  12.  (Rn,Rnd  Sequence  vs.  Time) 
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provides  a  comprehensive 


Over  this  time  interval,  the  sequence  Rp’ 

representation  of  the  field  scattered  from  the  target.  In  applying  the  spectral  estimation  algorithm  to  each 
data  block,  typically  one  assumes  an  order,  N,  of  the  number  of  scatterers  characterizing  the  target. 
Generally  one  chooses  N-large  and  then  filters  the  result  by  thresholding  on  the  amplitude,  |Dn|. 


Orders  1-5 


Figure  13.  Scattering  Amplitudes  (dB) 


In  Figure  12  we  illustrate  the  (Rp,  Rp)  sequence  assuming  N  =  10,  and  then  threshold  the  result 

for  the  five  dominant  centers  indicated  by  order  =  5  on  the  figures.  The  corresponding  |Dn|  (in  dB)  are 
indicated  in  Figure  13.  Noting  the  nearly  20-dB  reduction  in  lower  order  poles,  we  use  only  the  dominant 
five  poles  for  the  following  simulation. 

As  discussed  in  detail  in  Reference  1,  a  key  advantage  of  the  “snapshot  imaging'’  technique  is  that 
when  contrasted  to  other  Extended  Coherent  Processing  (ECP)  techniques,  only  a  limited  set  of  snapshots 
are  required  to  characterize  the  target.  Thus,  once  a  specific  scattering  center  is  identified  and  located  on 
the  target,  additional  data  does  not  add  significantly  to  the  location  estimate.  (Other  than  for  noisy  data, 
where  additional  data  increases  the  effective  SNR)  For  our  purposes,  to  illustrate  the  methodology 
developed  in  Section  6,  we  choose  a  subset  of  (Rp,  Rp)  data  over  the  time  interval  (0.43,  0.47)  seconds. 
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Over  this  interval,  the  target  aspect  is  nearly  constant  at  6  ~  60°,  |6|  «  |(j)|  and  the  spinning  motion 
of  the  target  dominates.  This  localized  motion,  in  effect,  looks  like  a  spinning  target  at  constant  aspect 
angle.  Figure  14  illustrates  the  (R^,  Rn)’  n=l,...,5,  sequence  over  this  interval  (vs.  sequence  label  q, 

where  t  =  ti,...tq.  and  q  ranges  from  (1  -  4533)  over  the  (0.43,  0.47)  time  interval  chosen.  We  also  plot 
the  (0,  (}))  variation  over  this  data  set  interval.  As  discussed,  9  is  nearly  constant  at  60°  and  (j)  ranges  over 
a  roughly  1 90°  interval  (115°-  305)  degrees. 


E 


- ► 

Time  Sequence  Label 


Figure  14.  (Rj^,  R^^and  G,  (j)  Used  for  Imaging  Simulation 


Figure  15  illustrates  the  results  of  applying  the  surface  constraint  imaging  algorithm  to  a  small  data 
subset  of  275  angle  steps,  over  which  ^  ranges  over  a  47°  sector,  and  this  sector  is  moved  progressively 
through  190°  to  cover  one-half  of  target  revolution.  The  scattering  centers  detected  are  plotted  in  "Ted"' 
overlaid  on  the  ‘"green’'  surface  model.  As  expected,  dominant  scattering  occurs  from  the  cylinder  edges 
and  the  fins.  Most  interestingly,  the  scattering  from  each  of  the  two  visible  fins  indicates  that  significant 
scattering  occurs  from  all  parts  of  the  fins.  Also  note  that  the  location  of  the  point  scatter  contributions  is 
consistent  with  incident  angle.  Figures  15b),  c)  and  d)  continues  the  sliding  47°  sector  image  sequence, 
progressively  showing  how  differing  scattering  centers  appear  on  the  target  as  the  target  rotates. 
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(t>  =  (113,160)  Deg 


<t)  =  (160^  207)  Deg 


<t)  =  (254^  30r)  Deg 


(|)  =  (207,  254*^)  Deg 


Figure  15.  Sequential  47 -deg  Sector  Images 


Finally,  we  illustrate  in  Figure  16  the  composite  of  the  images  in  Figure  15  -  i.e.  ^  ranges  over  the 
full  interval  (1 13°  -  301°).  We  show  a  3D  view  as  well  as  top  and  side  views.  Observe  that  we  do  pick  up 
a  few  stray  false  scattering  points.  These  however  result  from  the  output  of  the  spectral  estimation 
algorithm  (See  Figure  12)  as  sometimes  "‘glitch  points’"  occur.  We  have  yet  to  perfect  this  algorithm  to 
eliminate  these  points.  They  could  readily  be  filtered  out,  however,  using  a  "'noise  filter"’  in  series  with  the 
ambiguity  filters. 

Careful  examination  of  the  scattering  centers  extracted  from  Figure  16  indicates  that  a  component- 
based  signature  model  of  this  target  can  be  constructed  using  only  the  fm-components  and  three  slipping 
scatterers:  the  visible  truncated  cylinder  forward  and  rear  slipping  scatterers  and  a  slipping  scatterer  at  the 
visible  rear  base  of  the  cylinder.  The  coefficient  D„  (0,  (j))  for  each  of  these  components  is  readily 

determined  from  the  extracted  pole  complex  amplitudes  (see  Figure  13),  and  the  scattered  field  can  be 
reconstructed  by  superposition  of  field  scattered  from  each  of  these  scattering  components. 
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Top  View 


3D  View 


Figure  16.  Composite  Image  Points  Over  190°  Spin  Angle 
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8*  SUMMARY 


In  this  paper,  we  extend  the  novel  formulation  for  3D  imaging  of  ISAR  data  using  recent 
developments  in  high-resolution  spectral  estimation  theory  to  the  situation  where  some  apriori 
information  about  the  target  surface  geometry  is  known. 

For  some  applications  one  has  apriori  information  about  the  target  geometry  -  For  example  in 
processing  static  range  data  or  in  processing  flight-test  field  data  on  known  targets.  Of  prime  interest  in 
this  case  is  identification  of  the  dominant  scattering  centers  on  the  target,  with  the  intent  to  form  a 
computationally  efficient,  measurements  based,  signature  model  of  the  target.  In  this  report,  we  develop  a 
methodology  which  incorporates  this  known  surface  geometry  constraint  to  generate  a  3D  ISAR  image 
identifying  the  principal  scatterers  contributing  to  the  target  radar  cross-section.  The  methodology  uses 
only  a  sequence  of  single  snapshots  of  data  from  which  key  features  on  the  target  can  be  readily  extracted. 
The  technique  eliminates  the  need  for  a  range-Doppler  correlator  required  to  track  any  specific  scatterer. 
This  is  accomplished  by  hypothesizing  a  given  target  surface  constraint  (3D),  and  mapping  sequential  2D 
range-Doppler  location  estimates  onto  this  surface.  The  process  is  adaptive,  as  the  assumed  geometry  is 
modified  as  appropriate  to  fit  the  data.  When  applied  to  the  sequence  of  high-resolution  snapshots 
obtained  using  the  direct  spectral  estimate  obtained  from  sequential  blocks  of  data,  it  presents  a 
computationally  efficient  way  of  generating  a  high-resolution  3D  image  where  the  dominant  scattering 
centers  on  the  given  surface  constraint  are  identified  and  characterized. 

We  have  presented  several  examples,  which  illustrate  the  utility  of  the  technique.  Two  generic 
shapes,  a  cone  and  a  cylinder,  allow  for  an  analytic  solution  leading  to  a  quadratic  equation  characterizing 
the  scatterer  locations.  The  two  roots  obtained  from  the  quadratic  correspond  to  a  true  scatterer  location, 
and  an  ambiguous  solution  that  must  be  filtered  out.  We  consider  four  physical  mechanisms  for  these 
ambiguity  filters:  shadowing;  slipping  scatterers;  component-to-component  coupling  and  CAD  model 
surface  tolerances.  The  latter  two  are  applicable  to  the  case  considered  where  we  apply  a  more  general 
surface  constraint  characterized  by  Eq  (32).  We  illustrate  these  filters  applied  to  a  more  complex  target 
geometry  consisting  of  a  finite  cylinder  having  four  equally  spaced  fins  at  the  rear. 

Applications  of  the  technique  illustrate  the  snapshot  nature  of  the  component  extraction  -  i.e.  one 
not  need  a  continuous  wide-angle  sector  of  data  to  generate  a  3D  image.  Thus  it  is  particularly  applicable 
to  3D  imaging  of  sectors  of  sparse-angle  data. 
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The  extracted  scattering  locations  (Xn,  yn,  Zn)  and  component  features  (i.e.,  |D„|  vs.  angle)  can  be 
used  to  develop  a  computationally  efficient,  measurements-based  RCS  signature  prediction  model.  A  key 
attribute  of  the  resulting  model  is  that  the  measured  RCS  is  directly  incorporated  into  the  computational 
model.  The  extracted  scattering  locations  and  component  features  can  be  combined  using  the  Geometrical 
Theoiy  of  Diffraction  (GTD)  to  reconstruct  the  scattered  field  in  the  data  observation  region  (see  Eq  (2)), 
and  using  physics-based  characterization  of  the  components,  extrapolated  outside  the  observation  region 
for  wide-angle  characterization  of  the  target  RCS. 
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